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Abstract 



^ , 

^ ' Recently, it has been proved that a nonlinear quantum oscillator, gen- 

O^. eralization of the isotonic one, is exactly solvable for certain values of 

its parameters. Here we show that the Schrodinger equation for such an 

oscillator can be transformed into a confluent Heun equation. We give a 

^ ■ very simple and efficient algorithm to solve it numerically, no matter what 

f--^ ' the values of the parameters are. Algebraic quasi-polynomial solutions, 

^SJ , for particular values of the parameters, are found. 

(N 

x^ '. 1 Introduction 

O, 

f^ ^ Two years ago, Carinena, Perelomov, Ranada and Santander [Ij considered a 

quantum oscillator, intermediate between the harmonic and the isotonic ones, 
whose Schrodinger equation reads 



d 



a 






.'2J2 



w^x^ + 2g 



(x2 + a^f 



* + 2£;* = 0. (1) 



The interest of those authors on that problem lay on the fact that, as they 
proved, it is exactly solvable for certain values of the parameters, namely, 5 = 2, 
and w and a such that ujo^ = 1/2. Very recently. Fellows and Smith |2] have 
shown that this particular case of generalized isotonic oscillator is a super- 
symmetric partner of the harmonic oscillator. This fact has allowed them to 
reproduce the results of Ref. jl] in a very concise and elegant manner, and even 
to construct an infinite set of oscillators, with potentials approaching that of 
an isotonic oscillator as a; — >■ oo, all of them being partners of the harmonic 
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oscillator and, therefore, exactly solvable. Related also to the generalized iso- 
tonic oscillator of Carifiena et al. , another class of exactly solvable problems has 
been obtained by Kraenkel and Senthilvelan [3] . They have used point canonical 
transformations to convert Eq. ^ into a Schrodinger equation with a position 
dependent effective mass, which, with adequate mass distributions, may rep- 
resent different problems encountered in semiconductor physics. On the other 
hand, in view of the considerable progress in the synthesis of artificial quan- 
tized structures, a great variety of shapes for the potential wells and barriers 
are easily feasible and it seem plausible to see generalized isotonic oscillators as 
possible representations of realistic quantum dots. 

The rising interest on generalizad isotonic oscillators has lead us to try to 
contribute to a better understanding of their main features. Dealing with Eq. 
([T]) without any restriction, apart from the trivial ones 

w > 0, g > 0, a^ > 0, 

about the values of the parameters, we have found that a very natural change of 
variable transforms it into a confluent Heun equation (CHE). Then, the energies 
E can be easily obtained as the zeros of a function defined by a continued 
fraction. Of course, the results of Ref. [1] are reproduced by our procedure. 
This allows us to relate, by a continuous variation of the parameter g, the 
energy levels of the harmonic oscillator {g = 0) with those of the generalized 
isotonic oscillator. Besides this, quasi-polynomial (i. e. product of a rational 
function times an exponential times a polynomial) wavefunctions appear for 
specific energy levels and particular values of the parameters ui, g and a, different 
from those considered in Ref. [1]. 

In Section 2 we transform the Schrodinger equation ([T]) into a CHE. To 
solve it, we propose an extremely simple algorithm that allows us to obtain 
the eigenvalues and the eigenfunctions with a great accuracy. Results of that 
algorithm are shown in Section 3. quasi-polynomial solutions of the Schrodinger 
equation are obtained in Section 4. The relation of the polynomials found by 
Carifiena et al. [T] with confluent Heun polynomials is discussed in Section 5. 
Some final comments are added in Section 6. 



2 A confluent Heun equation 

Equation ^ presents two regular singular points, at a; = ±m, and an irregular 
one of s-rank 3 at infinity. (Along this paper we adopt the definition of s-rank of 
an irregular singular point given in Sec. 1.2 of the book by Slavyanov and Lay 
[3].) Taking advantage of the symmetry of the potential, the number and/or 
rank of the singularities can be reduced by a very natural mapping, namely, 

z = x^/a^ . (2) 

In this way, the upper half of the x-plane goes into the whole z-plane and the 
real axis x e (— oo,-|-oo), relevant from the physical point of view, is mapped 



into the positive real semiaxis in the z-plane, covered from cx) to along the ray 
argz = 27r and from to oo along the ray argz = 0. The Schrodinger equation 
turns into 



d2* 1 d* 

dz^ 2z dz 



r f ,2^4 


q z-1 Ea^ 
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2 z(z + l)2 2z J 



* = 0, (3) 



where the two free parameters tua^ and g are assumed to be given, whereas Ea? 
represents the eigenvalues to be determined. That equation has two regular 
singularities, at —1 and 0, and an irregular one of s-rank 2 at infinity. This 
suggests to compare it with some one of the various forms of the CHE, which 
presents the same pattern of singularities. Extensive discussions of the CHE 
can be found in Refs. [4| and [5 . Equation ([3]), when written in the form 

*-0, (4) 

is an example of CHE in its natural form [3 Eq. (1.1.4)]. The singularity at 
z — —\^ coming from the singularity at a; = ia, has indices 

1 / , -\ 1 



d^-^ 1 d^ 

dz^ ' 2z dz + 


■ Lu^a'^ 5/2 
4 z 


5/2 
z + 1 


9 Ea^] 


(z + l)2 2z J 



Pi = 2 (l + v/1 + 4gj , P2 - 2 (l - Vl + 4.9) • (5) 

The mapping ^ has introduced the singularity (branch point) at z = 0, with 
indices 

z/i =0, 1^2^ 1/2 , (6) 

which correspond, respectively, to even and odd solutions "^(x). The singularity 
for z — )• CX) corresponds to that for a; — > 00. It is immediate to check that the 
two Thome formal solutions of ^ for z — > cx) behave as 

exp (± ^ z^ z^± (1 + 0{z-^)) , with ^± = T^ - i . (7) 

The problem of finding the energy levels of the generalized isotonic oscillator 
reduces, thus, to determine the values Ea^ such that the Floquet solutions (series 
of increasing powers of z) of ([S]), corresponding to the indices i^i or 1^2, vanish at 
infinity according to ([7]) with the minus sign in the argument of the exponential. 
This is the well known connection problem of the singular points at and at 00. 
The presence of the singular point at z = — 1 prevents the convergence of the 
Floquet solutions out of the unit disc. So, a process of analytic continuation 
would be necessary to find their behaviour for z — > cxd. Instead of this, Slavyanov 
and Lay [U Sec. 3.6] suggest to carry out a Jaffe transformation consisting in an 
adequate linear transformation of the dependent variable followed by a Mobius 
transformation of the independent variable to convert the interval z e [0, 00) into 
the interval [0, 1] for the new variable. We have found convenient to substitute 

*(z) = (z + l)^ exp(- — zj ii;(z), with /i = /i_ = _ - - , (8) 



in © to get 




d'^w ( y 1 2n \ dw 
dz^ \ 2z z+lj dz 


+ 


^^i{Lua^ + l/2) ^l{^i-l) g z-l ' 
[ z{z + l) (z + l)2 2 z{z+iy\ 



w = 0, (9) 

and then to apply the Mobius transformation prcseribed by Slavyanov and Lay 

t = ^— , (10) 

z + 1 ' ^ ' 

which transforms Eq. Q into (keeping the same symbol to represent the de- 
pendent variable in terms of the new independent one) 



w f Loa^ 1 2/i — 2 \ dw 



dt^ V (1 - ^) 2i(l -t) l-t J dt 
Xwa2 + l/2) ^ii^J.-l) g 2t-l 



t(i-i)2 [i-ty 2t{i~ty 



w = 0, (11) 



to be solved in the interval t € [0, 1]. The energy levels are the values of E such 
that the "even" [u = vi = 0) or "odd" {v = V2 — 1/2) series solution 

oo 

u;(i)=r^c„t", co^O (12) 

n=0 

becomes finite at i = 1. Substitution of (fT2)) in (fTTj) gives for the coefficients c„ 
the recurrence relation 

^n+l C„+i + Bn C„ + C„_i C„_i — 0, (13) 

where we have abbreviated 

A„j ~ {m + v){m + V — 1/2) , 

Bm = [m + v){-2m - 2v - ujo^ +2^j.- 1/2) + n{ua^ + 1/2)+ g 12, 

Cm — {'m + v){ni + l + v-2^L)+^i{ji.-l)-g. 



The recurrence relation ([T3|) is an irregular difference equation of the Poincare- 
Perron type. The procedure suggested by Slavyanov and Lay to solve the con- 
nection problem consists in considering the physically acceptable solution of 
(fT3|) as a linear combination of its two Birkhoff solutions and adjust the value of 
Ea^ so as to reach the cancelation of the coefficient of the exponentially diver- 
gent one. Instead of this, we have preferred to base ourselves on the algorithms 
proposed by Gautschi [6] to find minimal solutions of three-term recurrence re- 
lations. The crucial consideration is that the solution minimal for n — ^ c» turns 
out to be dominant when the recurrence is used "from tail to head" . Bearing 



this in mind, we write the above recurrence relation in the form 

Cn — 1 J^n (-^n—1 



(14) 



^n ^^n — 1 

Cn+1 

and use it to compute c_i/co, starting with the approximate value 

— ^— ~ 1 H -J— , for n sufficiently large, (15) 

c„+i nV^ 

obtained from the characteristic equation |4j Sec. 1.6. 3] of the recurrence relation. 
The energy levels are obtained from the values of Ea^ such that one gets 

C-i/co = 0, (16) 

which implies c_i = and, in view of (|13|) . c_2 = c-3 = . . . = 0. 

3 Some results 

We have applied the above described algorithm to the determination of the 
four lowest states of a generalized harmonic oscillator with intensity g varying 
from g — (harmonic oscillator) to g — 20, for two typical values of wa^, 
namely tua'^ = 1/2 and ua'^ — 2. The results are shown in Figures 1 and 2. The 
behaviour of the eigenenergies, as g increases from zero, can be easily understood 
in view of the probability density of the harmonic oscillator eigenstates and the 
fact that the additional potential 



^ (a;2 + a2)2 

is negative for |a;| < a and positive for |a;| > a pj Fig. 1]. For low values of g, the 
effects on the energy of the positive and negative parts of that potential almost 
cancel to each other, except for the fundamental state whose density probility 
concentrates near the origin, where the additional potential is negative. This 
explains the gap between the energies of the fundamental and the first excited 
states encountered by Carinena et al. [1] in the case of ujo? = 1/2 and g = 2. 
As g increases further, the deeper and deeper potential well at the origin makes 
to decrease the energies of more and more excited states. 

For possible numerical comparison with results obtained by other methods, 
we report ours, for some arbitrarily chosen values of g, in Tables 1 and 2. The 
energies of the four lowest states, given with ten decimal digits, have been 
obtained by using a double precision FORTRAN code. The presence, in Table 
1, of exact values of the energy in the case of g = 2 is not surprising: this is the 
exactly solvable case discussed by Carifiena et al. ^Lj. More intriguing are the 
exact values of Eio^ and E2a^ found in Table 1 ioi g — 12 and those of Eoo^ 
for g = 20 and of Eia^ for g = 42 in Table 2. As we are going to show in the 
next Section, these are also cases of quasi-polynomial wave functions. For those 
values of the parameters, one could speak of quasi-exactly solvable potentials. 




Figure 1: Variation of the four lowest energy levels of the generalized isotonic 
oscillator with the intensity g. The parameters a; and a of the oscillator are 
assumed to be such that ojo^ — 1/2. 



Table 1: Energies of the four lowest states of the generalized isotonic oscillator 
of parameters lu and a such that cja^ ~ 1/2 and for some particular values of g. 
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Eoa^ 


Si a" 


£^20' 


Esa' 


0.00001 


0.2499965567 


0.7500003296 


1.2499993408 


1.7500000608 


0.1 


0.2145449837 


0.7531700367 


1.2438631419 


1.7505612654 


1 


-0.1849646064 


0.7686368018 


1.2241726194 


1.7519905648 


2 


-0.75 


0.75 


1.25 


1.75 


10 


-6.6939292571 


-1.3426308508 


1.4335481682 


2.1092577781 


12 


-8.3296211722 


-2.25 


1.25 


2.2548656278 


20 


-15.1074137799 


-6.5820180239 


-0.6097201160 


2.6043993900 


50 


-41.9395696487 


-27.0750778124 


-14.7363177369 


-4.9623193955 




Figure 2: Variation of the four lowest energy levels of the generalized isotonic 
oscillator with the intensity g, for cua^ = 2. 



Table 2: Energies of the four lowest states of the 
of parameters lu and a such that ujo^ = 2 and for 



generalized isotonic 
some chosen values 



oscillator 
of 5- 
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Eoa^ 


Eia' 


E2a' 


Esa' 


0.00001 


0.9999937095 


2.9999977428 


4.9999984646 


6.9999989879 


0.1 


0.9368657901 


2.9772742737 


4.9847133541 


6.9898929491 


1 


0.3495953307 


2.7588911779 


4.8519466428 


6.9003013951 


2 


-0.3372372644 


2.4870257918 


4.7099762556 


6.8039923347 


5 


-2.5490351910 


1.4941832183 


4.2680431727 


6.5346852493 


10 


-6.5291427792 


-0.6609393149 


3.3184939783 


6.1004000480 


12 


-8.1825461552 


-1.6592922308 


2.8380146272 


5.9058815492 


20 


-15 


-6.1825461552 


0.3407077692 


4.8380146272 


42 


-34.5910045651 


-21 


-9.8551888685 


-1.1281494657 


50 


-41.8769597362 


-26.8630723075 


-14.3102873432 


-4.2061920738 



4 Quasi-polynomial solutions 

The study of the generalized isotonic osciUator done by Cariiiena et al. [I] 
revealed the existence of quasi-polynomial solutions for certain values of the 
parameters. Specifically, they found that, if the parameters are related in the 
form 

g = 2uja^{2uja^ + 1), (17) 

equation ([T]) admits a quasi-polynomial solution 



*o - (^2+^2)2.^^ exp [-^^x'j , E, = -u:- i2^ar, (18) 

A'o being a normalization constant. It corresponds to the ground state of that 
generalized isotonic oscillator. More interestingly, they found that if, besides 
Eq. ((ni), one has 

1/2, (19) 



2 



all eigenstates, of energies 

Em = {m- 3/2) w, 
are represented by quasi-polynomial wave functions 

■^m{x)^ , ,"' ,. exp(-^x")7^„,(c^^/'x), m = 0,3,4,..., (20) 

a;(a;^ -I- a^j V 2 / 

the polynomials Vm being linear combinations of three consecutive Hermite 
polynomials of the same parity. 

The fact that simple analytic solutions of quantum problems, for particular 
values of the parameters, can serve as a check on numerical calculations was 
already pointed out by Demkov [7] in his study of the motion of a particle in 
the field of two Coulomb centers. He found that for certain set of values of the 
charges Zi and Z2 of the centers and the distance R between them, the wave 
function becomes quasipolynomial when written in ellipsoidal coordinates. 

In order to analyze the possible existence of quasi-polynomial solutions of 
Eq. ([3]), it is convenient to make explicit the behaviour of such solutions at the 
singular points. With this purpose we write 

vl/(z) = (z + l)''z''exp('-^z')$(z), (21) 

where p represents either pi or p2, given by ([5]) and such that 

P(P- 1) =3, 

and, according to ([S]), i' = J^i = in the case of even wavefunctions and v = 
V2 — 1/2 for odd ones. The equation satisfied by $(z), 

(P^ ( , 2z^ + l/2 2fl \ d$ 
az'' \ z z + 1 J dz 

ujg^ifi -V- p)z + oja^jp - ly) + p{2v + p/2) 
^ z{z + l) ' ^ ' 



with fi given in Eq. (|S]), is again an example of CHE. In fact, the change of 
variable z — > —z permits to compare it with the non-symmetrical canonical 
form of the CHE [SJ Eq. (1.2.27)]. Polynomial solutions of the CHE have 
been studied by Slavyanov Sec. 3.4]. His analysis, however, is not directly 
applicable to our problem due to the different role played by the parameters in 
his equation and in ours. 

It is possible to write a formal solution of (|22l) as a power series 



Hz) 



n=0 



with coefficients obeying the recurrence relation 

Ckn+2 an+2 + Pn+1 fln+l + "fn a-n = 0, 

with 



(23) 



(24) 



fl-i = 0, 
where we have abbreviated 

= m(m + 2j/-l/2) 



flo ^ 0, arbitrary. 



m,{m ~ uja^ + 2u + 2p - 1/2) + uja'^{fi - j/) + p{2i^ + p/2) 



reduces to a polynomial 



7„ = -uja^{ni- {p-v - p)) . 
Obviously, the series in the right hand side of 

k 

Hz) = Qk = ^a„,fcz", 

n=0 

the coefficients a„^fc being solution of (|24l) . if the two conditions 
7fe = and afc+i,fc = 



(25) 



(26) 



are satisfied. Assuming that the parameters uj and a? are given, the first one of 
those conditions, gives the eigenenergy 



E = {2k + 2p + 2v + l/2)L 



(27) 



in terms of p, whereas the second one determines the values of p, and conse- 
quently of g, for which the quasi-polynomial (PT|) . with $(z) replaced by Qk{z)^ 
is a solution. Notice that the second condition ([25)1 can be equivalently ex- 
pressed as the cancelation of the determinant of a tridiagonal (fc + 1) x (fc + 1) 
matrix. 



det 



/ /3o 


ai 
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70 


/3i 


a2 








71 


P2 


as 

7fc-2 


h~i ttfc 


v 








7fc-i Pk / 



(28) 



with fi replaced hy k + v + p iw the expressions of am, Pm and 7,„ given above. 
The left hand side of p8)) is a polynomial of degree 2{k + 1) in p. One of its 
roots, for every value of k,\s p — 0: the pure harmonic oscillator possesses quasi- 
polynomial solutions of every degree. Obviously, conjugate pairs of complex 
roots may appear, but these are not interesting, in view of the restriction g > 0. 
The real roots are easily obtained algebraically or numerically. To illustrate, we 
report the results obtained in the cases of fc = 0, 1, 2 and for the two values of 
Loa^ considered in Tables 1 and 2, respectively. 

4.1 Case k=0 

The first of conditions (PE)) . 70 = 0, gives p — v + p and, in view of (|H]), 

E ^ {2p + 2v + 1/2) uj. (29) 

The second condition, /3o = 0, gives 

p {ujo^ + 2v + p/2) = 0, (30) 

that, besides the trivial one p = 0, has the solution 

p = -(2wa^+4i/), i.e. g ^ {2uja^ + ^v){2uja^ + Av + I) . (31) 

4.1.1 Even solutions {v = 0). 

In the case of being 

p = —2uja , i. c. g = 2uja {2uja + 1), 

and for the energy 

E = {-Aloo^ + 1/2) uj, 

one has the quasi-polynomial solution 



*(z) o^{z + l)-^^"- exp ( — z 

or, in terms of the original notation, 

^{x) ex {x^ + a2)-2<^«' exp (--x^) ■ 

This is the solution mentioned in Eq. (2) of Ref. li . It can be recognized in our 
Tabic 1 (ground state for 5 = 2) and in our Table 2 (ground state for g ~ 20). 

4.1.2 Odd solutions {v = 1/2). 

For 

p = -2a;a^-2, i.e. g = {2uja^ + 2){2oja^ + ?,), 
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and energy 

E = (-4wa^-5/2)w, 

one finds the quasi-polynomial solution 

*(z) oc (z + l)-2-^-2 ^1/2 exp (-'^z 

or, in terms of the variable a:, 

vI/(x)cx(a:2+a2)-2--'-2a;cxp(-|a;2) . 

Examples of the case under consideration appear in Tables 1 (first excited state 
for g = 12) and 2 (first excited state for g = 42). 

4.2 Case k = 1 

4.2.1 Even solutions (v = 0). 

Now, from (P7|) . we have 

E^ (2p + 5/2) Lu. 

and, from (|28|) . 

- p (p^ + 4(wa2 + l)p^ + {Wo" + lOcjo^ + \)p + 2wa2(2wa2 + 5)) = , 
For ujo} = 1/2 there is a nontrivial real solution 

p = -^l/3/3_2-5^"^/^ with A = 3 fse - \/92TV 
corresponding to an intensity of the additional potential 

g = ^2/3/9 + 5^1/^3 + 28/3 + 25A"^/^ + 25^"^/^ . 
For wa^ = 2 there is also one nontrivial real solution 

p = -5^/^3-4- 11 B"^/^ with B = 3 ^72 - \/Ti9TV 

to which it corresponds 

g = 5^/3/9 + 35^/3 + 82/3 + 995^1/^ + 121B"2/3 _ 

4.2.2 Odd solutions {v = 1/2). 
For energy we have the value 

E^ {2p + 7/2) u. 
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and for p the equation 

- p (p3 + 4(cja2 + 2)p2 + {^uj'^a'^ + Igcja^ + 15)p + 2(2^^0'^ + Qua^ + 3)) = , 
For uja? = 1/2 one gets three nontrivial real solutions, namely, 

p = -l, p = -i(9-A/T7), p=-l(9 + \/T7' 
with corresponding intensities 

5 = 2, 5 = 29-5\/T7, g = 29 + 5VT7. 

For Loa? — 2 one obtains also three nontrivial real solutions 

p = -(1/3) (16-55^/3 •23ffn-l + SiVeV^jj , 

p = -(1/3) ne + 55^/3 sftni-iVs") f-i + siVeV ^ 

The corresponding values of g are obtained by using 

g = p{p-l). (32) 

4.3 Case k = 2 

4.3.1 Even solutions {v = 0). 

The energy is now given by 

E^{2p + 9/2)u;, 
where p is a solution of 

- p(p^ + 6(oja^ + 2)p^ + 3{4Lu^a^ + ISuja^ + 13)p^ 
8 V 

+ 4(2^^0^ + ISoj^a"^ + 39oja^ + 8)p^ + 2(12u;-''a^ + 82uj^a'^ + lOSojo^ + 3)p 
+ 4(4^^0^ + 28uj^a'^ + 27uja^)] = , 

For iua^ = 1/2, this equation has, besides the trivial solution, two complex and 
the real ones 



P = -1, 
P 



= -^f7 + Ci/2_f37-C + 64C-i/^^'^' 



2 



p = -^n + c^^^ + Osi-c + eAC'-^/^y j , 
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where we have abbreviated 

C = i ( 37 + 5^/3 ( fyge? - 6 V883749^ ^ ^ + fTQe? + 6 VSSST^) ^ ) ) " 

The corresponding intensities of the additional potential are given by Eq. ([5^ . 
For wa^ = 2 there are also three non trivial real solutions, 

p ~ -1.0813491830140753 , p ~ -7.22694217991115 , p ~ -11.525122115065377 , 

which correspond to 

g ~ 2.2506652386192836 , g ~ 59.45563545169008 , g ~ 144.3535618822344 . 

4.3.2 Odd solutions (j/ = 1/2). 

We have for the energy 

E^ (2p+ll/2) uj, 

and for p 

I p(p^ + 6{uja^ + 3)p'^ + 3{4:uj^a'^ + 26uja^ + 35)p^ 
8 V 

+ 4(2a;3a^ + 24a;2a'^ + Sluja^ + 59)p^ + 2{Uuj^a^ + 118a;2a'* + 276wa^ + 105)p 
+ 4(4a;^a^ + iOuj'^a* + Ihuja^ + 15)] = . 

For ujq} = 1/2 we obtain three nontrivial real solutions, 
V = -1, 

1 / / / / \ 1/2 



where 



Z) = i U6 + (56575 - 12 VM56U7) + (56575+12^3456147) 



1/3 



The corresponding values of g follow from Eq. 

For wa^ = 2 there are also three nontrivial real solutions, 

p ~ -1.0527990898606965 , p ~ -9.18056490570385 , p ~ -13.30773725973479 , 

corresponding to intensities 

g ~ 2.1611850134722075 , g ~ 93.46333689354499 , g ~ 190.40360823493683 . 
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None of the above considered cases explains the value £^2^^ = 5/4 found for 
ua'^ = 1/2 (Table 1) and g = 12. It is not difhcult to see that it corresponds to 
fc = 4, 1^ = 0, and p = — 3. The wave function is in this case 

*(z) ex (z + l)-^ cxp(-z/4) (1 - 10 z - 4 z^ _ (2/3)z3 - (l/21)z'^) , 
that is, in terms of x, 

10a;2 4x4 2x^ x^ 



*(x) oc {x^ + a^)--^ exp(-a;74a-^) 1 - 



3a6 21aS 



5 The exactly solvable case 

In the preceding Section we have found, whenever wa^ = 1/2, a quasi-polynomial 
solution with p — —1 for almost every considered value of k and in both cases 
of even or odd wave functions. There are only two exceptions, namely the case 
oi k = 0, f = 1/2 and that oi k — 1, v = 0. Although we have not discussed 
the cases of fc ~ 3,4, . . ., one can easily check that the same value p — —1 
appears for z/ = and for f = 1/2 in all cases. This was to be expected, since 
for uia^ = 1/2 and p = —1, that is, 5 = 2, the generalized isotonic potential is 
exactly solvable [T]. In what follows, we show that the polynomials entering the 
solutions found by Cariiiena et al. are in fact confluent Heun polynomials. 

We assume from now on that ioa? = 1/2, p — —1, and p — k-\-i> — 1, k being 
a positive integer. According to Eqs. ((2T|) . ((22)) and (l25t . the quasi-polynomial 
solutions are of the form 

*(z) = (z + l)-iz'^exp(-z/4)Qi'')(z), i^ = 0, 1/2 , (33) 

the polynomial (of degree k) Qj^ obeying the differential equation 

d'Qi'\f_l 21^+1/2 _ 2 \ dQJ^^ fcz/2 + fc/2-2^ (.)_ 
~^^^[ 2+ ; VTl) ^T^ ^i^TT) ^' -U. (34) 

This is but a particular case of Eq. (1^^ that, as we have already mentioned, is 
a CHE. Confluent Heun polynomials can be written as linear combinations of 
hypergeometric or confluent hypergeometric polynomials 5, Sec. 2.3]. This sec- 
ond possibility is more convenient for a comparison with the results of Carihena 
et al. [1]. With this purpose, we introduce a new variable 

,-l_-^. (35) 

in terms of which the differential equation reads 

{y+l/2)y — ^ + {{y + l/2){2u+l/2-y) 2y) -^ 
ay ay 

+ {{y + l/i)k-2u)Q!^^^ =0, (36) 
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that can be written in the form 



((y + l/2)I?o + 2?i) Qi'^^^O (37) 

with the differential operators 

d^ d 

Do = y j^ + {2iy+ 1/2- y)— + k, (38) 

Vi = -2y^-2iy. (39) 

Now we try in ([37]) the sum of confluent hypergeometric polynomials 

fe 

Qk^ - E <fe ^^(-(^ - ")' 2^ + 1/2, y) , (40) 

ra=0 

with coefficients A^^l to be determined. According to Eqs. f3.f .f and f3.4.f0, 
respectively, of Ref [8) , one has 

VoM{-{k-n),2i' + l/2,y) = nM{-{k - n),2iy + 1/2, y) , (41) 

ViM{-{k-n),2i' + l/2,y) = -2{k - n + v) M{-{k - n),2v + 1/2, y) 

+ 2(fc - n) M{-{k -n) + l,2v + 1/2, y) , (42) 

and Eq. (|37)) turns into 

fe 



E ^S ( ("(2^ + 5/2) - 2fc - 21^) A^(-(fc - "): 2;^ + 1/2, y) 

+2(fc - n)M{-{k -n) + l,2v + 1/2, y)) = . (43) 

Cancelation of the coefficients of the successively decreasing powers of y in the 
left hand side of the last equation gives, for fc > 2, 

AqI arbitrary , 



^^^ " fc + 21. - 1/2 -^"''^ ' 



(u) _ 2^ - 3fc + Akv - 2v (^) 

2,fe " r,iu I o,, r/omrT"o7; Q7o^ -^o,fe ' 



2(fc + 2v- l/2)(fc + 2z/ - 3/2) ' 

3,fe ^ ■^4,fc 



.4^''^ - A^"'^ - - 



in both cases of ;. = or ;. = 1/2. Substitution of these expressions in (|40p 
gives 

QT = <!■ (^M{~k, 1/2, y) - ^Af (-fc + 1, 1/2, y) + J^M{-k+2, 1/2, y)^ , 

(44) 
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Qk^'^^AT (M{-k,3/2,y)-2M{-k+l,3/2,y) + ^-^M{~k+2,3/2,y)) 

(45) 
Now we can use the relations between confluent hypergeometric and Hermite 
polynomials [SJ Eqs. 22.5.56 and 22.5.57] H Sec. 10.13, Eqs. (17) and (18)] 

777 ' 

M{-m, 1/2, y) = (-1)™7:^^ H,m (V^) , (46) 

[2m)\ 



2(2m+l)! 



V^A/(-m,3/2,y) = (-l)" ^,^^', ,., g2,n+i (v^) , (47) 



and choose 



to get 

Qf = ^2fe(Vy) + skH2k-2{Vy) + ^H^k ~ 3)H2k-i{Vy) , (49) 

^1/2 ga/2) ^ ^^^^^ (^) ^ 4(2fc + l)H2k-i{Vy) + 4(2fc + l)(2fc - 2)H2k-3{Vy) ' 

(50) 
which are, respectively, the polynomials V2k and V2k+i of Ref. [IJ, the variable 
being 

V2;= ^=- = V^a;. (51) 

V2a 



In the above discussion of Eq. (|43)) we have left aside the cases of fc = and 
k = 1. For fc = and i^ = 0, Eq. (|43l) is satisfied irrespective of the value of 

^g Q, whereas for v = 1/2 its fulfilment requires ^q q = 0, which implies that 

the resulting Qq is identically equal to zero. For fc = 1, instead, Eq. (|43)) is 

satisfied, for v = 0, only if ^q { = All = and consequently Q{ identically 

equal to zero, whereas for v = 1/2 it only requires A[ i = —2Aq { . All this 
is in accordance with the two exceptional cases offc = 0, i^ = l/2, and fc = 1, 
t" = 0, encountered in Ref. 1 and mentioned at the beginning of this Section. 

6 Final comments 

Most of the quantum mechanical problems which admit a simple solution are of 
the hypergeometric class: the corresponding Schrodinger equation turns, by an 
adequate transformation, into a hypergeometric or a confiuent hypergeometric 
one. In this paper we have shown that the generalized quantum isotonic oscil- 
lator belongs to the Heun class. Other examples of physical problems of this 
class can be found in Chapter 4 of Ref. [4j . 

Up to now, the most appealing feature of the generalized quantum isotonic 
potential was its exact solvability for certain values of its parameters |I] . The 
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fact, shown in this paper, that it is quasi-exactly solvable when those parameters 
take particular values spread over a wide range makes it to be especially suited 
to serve as a workbench to test the accuracy of approximate (perturbative, 
variational, etc.) methods of solution of the Schrodinger equation. 

Quasi-exact solvability of quantum Hamiltonians is closely related to their 
Lie-algebraic properties [TO]. A discussion of this topic with reference to the 
generalized isotonic oscillator would be necessary, but it lies out of the scope of 
this paper. 

The sequence {Vn} arising in the exactly solvable case [Ij does not include 
a linear (n = 1) nor a quadratic (n = 2) polynomials. Therefore, it cannot be 
used as a basis for an expansion. Similar sequences of polynomial eigenfunctions 
of a Sturm-Liouville problem have been found by Gomez-UUate, Kamran and 
Milson [n] and by Quesne [12]. In the case of the first authors, the sequences 
does not include the constant (n = 0) polynomial. Quesne has found sequences 
without the constant polynomial and also sequences without the constant and 
the linear polynomials. Nevertheless, the first authors have proved that such 
sequences, which they denominate exceptional polynomial systems, are a basis 
in their corresponding L^ Hilbert spaces. Besides, Gomez-Ullate, Kamran and 
Milson [TO] have given an extension of the Bochner's theorem applicable to those 
sequences of orthogonal polynomials, solution of a Sturm-Liouville problem, that 
start with a polynomial of degree one. It would be interesting to explore the 
possibility of an analogous extension of the theorem for sequences of polynomials 
like that discussed in Ref. [L, where the absent polynomials are not the lowest 
order ones. 
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